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ABSTRACT 

We propose a new method of producing neutrino pair beam that consists of a mixture of neutrinos and 
anti-neutrinos of all flavors. The idea is based on a coherent neutrino pair emission from excited ions in 
circular motion. High energy gamma ray much beyond the keV range may also be produced by a different 
choice of excited level. 
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1 Introduction 


Synchrotron radiation is a very useful tool of photon emission up to the X-ray energy range, providing a well 
collimated beam. We examine a similar problem of neutrino pair emission under a circular motion of ions. 
When excited ions with a high coherence are circulated, emission rates become large with neutrino energies 
extending to much beyond the keV region in the form of well collimated beam. Produced neutrino beam 
is a mixture of all pairs of neutrinos, including This gives a CP-even neutrino beam, hopefully 

providing an ideal setting to test fundamental symmetries of particle physics [ 2 ], in particular, to measure 
CP violating (CPV) phases in the neutrino sector [3], [2], [3], [5]. Circulation of highly stripped heavy ions 
is desirable to achieve the highest neutrino energy in the GeV region with the largest production rates. 

Our method of calculation may be adapted to synchrotron radiation that occurs at electron machine, 
giving essentially the same results as in [ 6 ], although our method of calculation is different. We shall make 
it clear how a GeV range intense beam of neutrino pairs is made possible if one uses excited ions instead of 
ions in the ground state. 

One may also produce high energy gamma ray much beyond the keV range by an appropriate choice 
of excited level of different parity. This may be very useful since the usual electron synchrotron can only 
produce the keV range photon. 

The rest of this paper is organized as follows. In the first two sections we shall explain our semi-classical 
approximation to treat the ionic motion as given classically and to calculate the probability and its rate of 
neutrino pair emission in the standard electroweak theory. In Section 4 we give the core calculation of a 
phase integral that appears in the rate calculation. We find that with excited ions the phase integral over 
time contains stationary points of the phase, leading to large neutrino pair emission rates. In the following 
section we compute the differential energy spectrum of neutrino pair production at synchrotron site. In 
Section 6 we discuss a similar problem of photon emission from electric dipole allowed atomic transition. 
When a good coherence among ions in the excited and the ground levels is prepared and maintained, it 
might even be possible to have a coherent gamma ray emission much like laser in the optical region. 

In a sequel paper we shall discuss neutrino oscillation experiments that can be done away from the 
synchrotron. 

Throughout this work we use the natural unit of ^ = c = 1. 


2 Semi-classical approximation 

The total wave function of a composite ion consists of a direct product of the central motion (CM) part of ion 
as a whole and its internal part as a consequence of separation of hamiltonian operator into an independent 
sum of two terms. For the neutrino pair emission process of internal atomic transition, |e) —>■ l^), the other 
CM hamiltonian part never contributes simultaneously. It only contributes to the cases of |a) —>■ |o), a = e, ( 7 , 
and this gives rise to the usual synchrotron emission in much the same way as in the electron machine. 
For the internal transition, the CM part of wave function gives a weight factor of its probability density 
(y the quantization volume) in the internal part of hamiltonian, in accordance to the general 
rule of the correct property of the lifetime under the Lorentz transformation oc 7 [7]. Here 7 is the boost 
factor of excited ion related to the constant velocity v of circular motion by u = — 1 / 7 ^ ~ 1 - 1 /( 27 ^). 

Strictly, one needs the instantaneous boost factor 7 (t) of a time dependent fnnction, bnt the emission region 
around the circnlar orbit is short, and one may replace this by the constant circular velocity. For the ion 
internal state we shall confine ourselves to a two-level system as an approximation, its ionic states being 
|e) and \g). The metastable state |e) in an upper energy level is assnmed to have the same parity as that 
of the ground state \g) such that fast electric dipole transition is forbidden, while a magnetic dipole (Ml) 
transition and the neutrino pair emission are both allowed. Relevance of the Ml transition to neutrino pair 
emission is explained in dne conrse. Another electric dipole case between different parity states is usefnl for 
high energy gamma ray emission and is discussed in Section 6 . 
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Bilinear forms of wave functions such as the ion current may be described in terms of the density matrix, 
Pab ,CL,b = e,g for its internal part. We assume that the central motion is described by the classical trajectory 
function x^it) (t the time at observation) of circular motion. The density matrix for two-level system is 
governed by the optical Bloch equation. Its solution may readily be derived in terms of initial values. In 
particular, the off-diagonal element Peg{t) of our system may be described to a good approximation [ 8 ] by 


Peg{t) = Peg(0)exp[-(ieeg 


1 , 

+ ’ 
^2 7 


( 1 ) 


when effects of photon emission are highly suppressed. Throughout this work we use the time in the 
laboratory system in which measurements of neutrino beam experiments are done. The phase relaxation 
rate l/r 2 is usually larger than its minimum value l/(2re) (re being the natural lifetime of state |e)) that 
occurs when the phase relaxation is dominated by the spontaneous decay. 

In the atomic physics community the quantity p^g is called the coherence. For a pure quantum state of 
a single atom, it is given by a quantum mixture of two states, |e) and \g). Its value is bounded to be less 
than the value 1/2 in our normalization convention. Its macroscopic average over a collective body of atoms 
or ions is usually much less than this maximum value. We shall not discuss the experimental problem of 
how a large initial coherence given by Pegi^) may be prepared. 

We can neglect contributions of ionic states that remain in either the excited or the ground state, their 
rates being proportional to p^^ , p^g, since they give rise to neutrino pair emission of much smaller rates and 
much smaller energies, the neutrino-pair analogue of the usual synchrotron radiation. This result originates 
from that these density matrix elements have no oscillating phase factor as in Peg{t) oc It is found 

that both in this case and in the case of electron synchrotron no large neutrino pair production occurs, as 
is made evident below. 


3 Perturbation theory of neutrino pair emission 

In our semi-classical approximation the hamiltonian system of interacting neutrino with atomic electrons 
is quadratic in neutrino field variables, and one can readily solve the problem of neutrino pair emission, 
using the perturbation theory of the weak coupling Gp- The four-Fermi interaction of neutrinos and atomic 
electrons is given by the hamiltonian (written in terms of neutrino field operators): 

+ cfjA>^{x))ul{x)apUj{x), (2) 

d ^2 jj = i 2,3 

with (ap) = (1, —(?). We use the neutrino index convention of Roman alphabets, a, b, c, to indicate neutrino 
flavor states, and Roman alphabets, to indicate mass eigenstates z^i,r' 2 ji^ 3 - The neutrino 

mass ordering is taken as usual: m 3 > m 2 > mi for the normal hierarchy case and m 2 > mi > m 3 for the 
inverted hierarchy case. Both W- and Z-boson exchange contributions are added, and the hamiltonian is 
written in the Fierz-transformed form (charge retention ordered). There are both vector and axial-vector 
currents, V{x),A{x), with their couplings ■ We may assume the non-relativistic limit for transitions of 
internal electron states in the rest frame of ion, which singles out as the dominant contribution the spatial 
part of axial 4-vector in the form of the electron spin current: (0 , Se ), Se = Yliai9\da/‘A^) 1 where the 

sum is taken over valence electrons of ions. Note that the monopole term of the vector part oc cTVq vanishes 
due to the orthogonality of wave functions between |e) and I 5 ). The coefficients of axial-vector parts are 

4 = c-c-t = l, (3) 

in the standard electroweak theory. 
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In the laboratory frame the relevant current becomes [9] 


(5„) = (7/3-5e,5e + 


7 


7 + 1 


(/3-5e)/?)~7(/3-^e,(/3-5e)/3) 


(4) 


where j3 is the Lorentz boost vector. Averaging over atomic spin direction to the leading order of large 7 
gives the squared amplitude summed over neutrino helicities m. 


si 


7'^ (1 + 7 


I Pi ■p2 mim2 


3 E 1 E 2 ‘^EiE2 


<5m 


(5) 


where 6m = 1 for the Majorana neutrino and 6m = 0 for the Dirac neutrino. Our experience of calculations 
for heavy atoms such as Xe, Yb, etc m, m shows that these matrix elements Se are of order unity or 
0(0.1) where the intermediate coupling scheme of heavy atoms holds. We assume in the following that the 
Ml transition matrix element is of this order. 

From this consideration it is found that the circulating atomic spin is the source current of neutrino pair 
emission and the relevant current is given by 


Jegi^) = S'^-^ J dtpegit) (x - XA{t) ) , 


= 


d^x^J^gix)- Cijv\{x)apUj{x), 


( 6 ) 

(7) 


where XA{t) = {t,fA{t )) (written in terms of the time in the laboratory frame) is the trajectory function of 
excited ion in circular motion given by 


^ , vt vt 

rAXt) = p sm — , 1 — cos — , 0 

p p 


( 8 ) 


where p is the radius of the circular orbit. The factor 1/^/7 in eq.dl]) arises from the overlap of CM wave 
functions, f d?X\^i{X)\^ = I/ 7 . 

We adopt the interaction picture in which the kinetic and the mass terms of neutrinos are taken as the 
free part of hamiltonian Hq, consisting of diagonal terms of did where b, d are annihilation operators of 
neutrino and anti-neutrino (in the Majorana neutrino case d = b) when they are mode-decomposed using 
plane waves of dehnite helicities. The Fermi interaction ([7]) due to the circular ion motion gives rise to 
off-diagonal terms, in particular, terms of the form, bd,b^dl In the perturbative picture this means that 
neutrino pairs may be created at ion synchrotron. 

The amplitude Aij{pihi,p 2 h 2 ]t) of neutrino-pair production of momentum p) (its energy given by Ei = 
\JVi S- rrif for neutrino of mass m,) and helicity hi time-evolves according to 


idtAij{pihi,p2h2;t) = idt{0\di{p2h2-,t)bj{pihi-,t)\0) = {0\[di{p2h2;t)bj{pihi;t), Hu,]\0), 


(9) 


reducing the calculation to the commutator between the neutrino bilinear field db and the weak hamiltonian 
Hyj. The result for neutrino pair emission of a single flavor is given by a time integral. 


1 If- 

Aij{pihi,p2h2;t) = -iy/2GF—Cij / j\{pi + p2]t') ■ ji^, 

yT J —00 

t) = Pegit)S°‘e~"d^'^^d) ^ = A{pihi)av{p2h2) ■ 


( 10 ) 

( 11 ) 


Here u, v are associated plane-wave solutions of emitted neutrinos. 

The basic interaction hamiltonian ([7]) indicates a number of striking features of neutrino pair emission 
process. Notably, it predicts a coherent (namely endowed with a dehnite phase relation among two neutrinos 
in the pair) mixture of all neutrinos and anti-neutrinos of three havors. 
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The semi-classical approximation in the present work is limited to the neutrino energy region in which 
ion recoil may be ignored, which allows the GeV region neutrino production since the circulating ion energy 
is much larger. 


4 Phase integral 

Pair emission rate defined by Pij{t]pihi,p 2 h 2 ) = dt\Aij{pihi,p 2 h 2 ',t)\'^ is given by 

Pij{to',Pihi,P2h2) = 4:Gp\peg{0)f-\Cij\‘^ [ dtS^^ {■N'ap{pihi,p2h2)P^^^^'^~^^^'> A 

= 4G|.|pes(0)p-|a,f r dtS^^ , (12) 

7 JO ^ ' 

A{t) = {El +E 2 - —)t - {pi + P 2 ) ■ rA{t) , (13) 

7 

AA"^(Pi/ii,P2 ^2) = juiPlhl,P2h2){jlf{Plhl,P2h2) , (14) 

by taking an infinite time limit, which effective means that time contributing to the integral (I12|) is much 
larger than a small fraction of the orbital period p/c. Explicit forms of Map{pihi,p 2 h 2 ) may be evaluated 
by using formulas given in m- 

The important phase factor in the integral is given by 

A(0) - A(-t) = {El + E 2 - — + i^)t - p({pi+ P 2 )x sin — + {pi+ P 2 )y (1 - cos — )] . (15) 

7 7^2 V P P J 

Let us introduce directional angles of emitted neutrinos: 

vt vt TT IT 

Pi =Pj(cosV’iCos(0j-I-),cosV’jsin(6'i H-),sin'0i), -tt < Oi < tt . (16) 

p p 2 2 

Angles are measured at an observation point away from the circular motion. The forward and the background 
directions with respect to the ion beam correspond to |0j| < tt/2 and |0j| > 7r/2, respectively. See Fig([T]) for 
this coordinate system. 

The phase factor A(0) — A(—t) of ea. (11511 contains three terms: in addition to the main term oc Ei + E 2 , 
one is from the circulating ion proportional to p, the radius of the orbit, and the other is proportional to 
the level spacing Cgg. Under the normal condition one may ignore the imaginary component oc I/T 2 , since 
€eg ^ 1/T2- The most important observation in the present work is that an input of de-excitation energy 
€eg may lead to cancellation of three terms and to existence of stationary phase points in the relevant phase 
integral along the real axis of time. On the other hand, without the €eg term one can show that the phase 
is positive definite. As is well known in mathematical physics, contribution around stationary points does 
not suffer from large suppression unlike constructive interference contributing with the same sign phase. 
This was the case without the €eg term such as synchrotron radiation and neutrino pair emission from the 
ground state ion. The well-known exponential cutoff arising from the constructive interference gives rise to 
the cutoff energy of emitted photon ~ 7^//o in synchrotron radiation [^. This cutoff also occurs for neutrino 
pair emission at electron synchrotron, restricting available neutrino energies up to a keV range. 

The crucial condition for the presence of stationary points is derived by setting vanishing time derivative 
of eq. (fT^ , leading to an equality, 

El + E 2 - — - Pi cos cos{6i + —) = 0. (17) 

^ i P 

Infinitely many stationary points exist along the real axis of time t. It turns out that the most important 
contribution comes from the point nearest to t = 0, the end point of the integration range t > 0. Note that 
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Figure 1 : Coordinate system for calculation of the phase integral. Observation is made at a point far away 
on the (positive side of) x-axis. Circular motion of excited ion in the (x, y) plane is depicted in red. The 
angle is defined to be zero in the ion orbit plane, while 6i is the angle measured from the tangential 
direction to the ion beam, with the negative region < 0 being defined towards the inner region of the 
circular orbit. 


without the e^g term there is no stationary point solution on the real axis, hence the possibility of large 
neutrino pair production at electron synchrotron is excluded. 

Let us consider the in-plane forward direction of V’i = = 0 for Ei S> €eg ■ One may expand the left hand 

side function of eq. (| 17 jl in powers of time variable t, to solve the stationary point condition. The look-back 
time of the stationary point t = tc is then derived as 


2e, 


eg 


7(.Ei + E 2 ) 7^ ’ 


(18) 


for El + E2 < 2 eegj. In the following we shall take the kinematical region in which this stationary point 
exists. In order to evaluate the phase integral to a good approximation, we use the power series expansion, 
to derive 


X = 


1^1 1 y 1/2 Pi COS ijji sin 6*1 -F P2 cos -02 sin 62 \ 

/2p V 7^ Pi cos V'l cos 61 + P2 cos "02 cos 62 ^ J ’ 


V2p 

-/I 

2 V 2 p 

l^r = (l 


A(0) — A(—t) ~ ^ ( —X — x^ 




- El - E2 + {1 — cos V’l cos 01 -I-P 2 COS ^^2 cos 02) + , 


2 = (1 - 1 ) 1/2 

T 


Pi cos V'l cos 01 -|- P2 cos V’2 cos 02 


^ — El — E2 + {I — cos 'tpl cos 01 -|- P2 cos '02 cos O2) + B 


1 1 1/2 (pi cos 01 sin 01-hp2 cos 02 sin 02)^ 

2 ^ 7^ Pi cos 01 cos 01 -|- P2 cos 02 cos 02 


(19) 

( 20 ) 
( 21 ) 
( 22 ) 
(23) 


In deriving this equation, we shifted the integration variable t such that 0 {t^) terms are eliminated. The 
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power series expansion in terms of time t has been retained up to 3rd order of because still higher order 
terms are suppressed by powers of tc/p- 

The condition that the stationary point is within the integration range t > 0 gives a limitation of emitted 
angles. An angular region deep inside the circle of ionic motion gives stationary points in the forbidden 
region of f < 0, hence does not give large neutrino pair emission rates. This forbidden region is defined by 

V2(.Ei0i + ^ 202 ) < + E2) - liEiiej + i/jj) + E2{9j + )(Si + E2) . (24) 

The necessary phase integral of x—variable involves a smoothly varying function h of time and it has a 
form, 

dx/i(x)cos^(^x^ - ~(Ji/3(0 + ^- 1 / 3 ( 0 ) > (25) 

where Juiz) is the Bessel function, and h(x) is a smoothly varying function of time given by squared matrix 
element of neutrino pair emission. The large radius limit of p ^ 00 , hence the ^(oc p) —>■ 00 , is important for 
calculation of differential rates, since the radius p is much larger than any microscopic length scale involved. 
The limit gives 



(ix/i(x)cosC(-x --x) ^ Y YCOs(^ --)-^ , (26) 

as ^ 00 . This asymptotic formula may also be derived directly using the steepest descent, or the 

stationary phase method of mathematical physics. The stationary point appears at x = 1, which implies 
that tc = \/2p/|A|. In addition to the phase factor given here there is a constant phase factor arising from 
the phase at the stationary point, namely A( 0 ) — A{—tc), which however gives a negligible contribution. 

The fact that the phase integral proportional to rate, ([2^ or (l26ll . can give negative values for some 
value of ^ might appear odd. But since this is time derivative of a positive quantity (probability), this 
may occur without any violation of fundamental principles. Indeed, this also occurs in the usual formula 
of synchrotron radiation [ 6 ]. The quantity ^ is actually a complicated function of neutrino energies, their 
emission angles, the boost factor 7 , and the atomic energy scale €eg- The region of these variables that 
effectively contributes with a large rate is found to give mostly positive rates. Thus, there is no serious 
problem of the negative rate. An alternative method used in the case of synchrotron radiation [ 6 ] treats 
the time and one of the angular variables, 9 (essentially not measurable) symmetrically in integration by 
changing integration variables in a clever way. A generalization of this method to the case of neutrino pair 
emission might be possible with much effort. 

For a finite value of T 2 , the stationary point moves to a point slightly off the real axis, introducing a 
small correction to €eg replaced by Cgg — The effect of this shift is small. 

It would be instructive, before proceeding, to mention the limiting case of €eg —>• 0 in our problem. In 
the limit the stationary point tc approaches the end point of time integration range and the phase space of 
neutrino momenta shrinks to zero. Thus, rate due to the mechanism considered vanishes in the limit. 

As another extension we would like to mention other contributions than the spin current contribution 
considered here. Contributions from the excited state and the ground state are proportional to pee, Pgg 
which do not have the Ceg factor in the phase, hence this case too has no stationary point on the real time 
axis. Result of the phase integral in these cases is given in terms of the modified Bessel function much like in 
the synchrotron radiation. The neutrino energy spectrum then suffers from the exponential cutoff of order 
7 ^/p, which is typically in the keV region. Since the weak interaction rates scales with energy^, rates are 
negligibly small. We have neglected these contributions. 
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5 Differential emission rate of a single pair 

Let us first write down the squared spin amplitude ([!]) in our coordinate system; 


1 


1 


1 


M. = - 7 ^<S'g 1 H— cos V’l cos 11)2 cos( 0 i — 62 ) -\— sin ijji sin 'tp 2 


( 27 ) 


For simplicity we took neutrinos to be massless which is adequate for our purpose. Effects of finite neutrino 
masses are significant only for Ei < O^mj), which is of order keV for m = 0 . 1 eV ,7 = 10^. Since neutrino 
pair production rates are small for this energy range, we shall ignore effects of finite neutrino masses for 
discussion of production rates. 

A straightforward calculation using this result gives the differential production rate for a neutrino pair 
I'ii'j of mass eigenstates. In order to avoid complication, we shall write down this formula in the leading 
approximation of large boost factor; 

(28) 


dEidE2dElidEl2 2’^/^ • 3\/^(27r)® 

. Epn El + E 2 . 1 

~ 2 


F={Ei + E2)C-^ - + Eli^l) - ^{9i - 92? - ^(Ei02 + E 29 I). (29) 


7 


^7 


The spin factor is given by A4 ~ in this approximation. The function F is more complicated in the 

most general case of the boost factor, which may be inferred from ea. (l44l) for (the electric dipole) photon 
emission. There is a constraint on angle factors given by F > 0. This constraint gives angular restriction 
worked out for apertures (given for simplicity to the case Ei = E 2 = E), 


= O ( - 


1 2{Em-2E) 


E 


A9 = 0 


Em ~ 2E 


En 


Em = 276 


eg ■ 


While A9i is of order unity individually, the opening angle of two neutrinos of the pair is limited by 

A\9i-92\ <0 


1 / {El + E2){Em — El — E2) 


7 


E 1 E 2 


(30) 


(31) 


The suppression by I /7 for the opening angle is of great interest from the point of oscillation experiments, 
since it leaves open for the possibility of a coherent neutrino pair interaction at measurement sites. We shall 
have much to discuss in the next section. 

Integration over four angle factors, with dflj = cos ipid'ipid9i, may be carried out, to give 


J dQi j ~ 


(Fi+F2)5/4^^ Fi+F2 ^7/4 

Em. 


(FiF2)3/2 TP > ^ (22) 

where V 4 = 7 r ^/2 is the volume of 4-dimensional sphere of unit radius. Using this result, the double 
differential energy spectrum becomes 


d?T. 




1 




dEidE 2 21 • 27 • 23/4^3^' 


TrTT’ 


7 


Er, 


(33) 

The relation = 1/4 was used. Further integration gives the single neutrino energy spectrum and 

the total pair production rate: 


dE 


1 


1 E 


21 • 240 . 


vrTT^ 


-StN\p,g{Q)\^^gG^pEl,-f{ 


7 E„ 


rl—x rl 

f{x) = i/x / dyy^^‘^{y + x?^*{l — X — y?^^ , / cix/(x) ~ 0.00727, 

Jo Jo 


0.0073 


21 • 210 • v^7r47 


S!N\p,g{ 0 )\‘^^GlE^„ 


(34) 

(35) 

(36) 





















Production rate does not depend on the mass eigenstate label, since we assumed the massless neutrino for 
this calculation. The normalized universal spectrum function f{x)/ dyf{y) is plotted in Fig([2|). The 
average neutrino energy is ~ 0.30Em- The end point behavior at a: = 1 gives the threshold behavior 
oc {Em — Ey^/‘^, Em = 2 eeg 7 the highest neutrino energy and oc y/E in the infrared region of PI —)• 0 . 



Figure 2: Normalized universal spectrum function f{x) of ea. (l35]) divided by its total integral dyf{y) in 
which X = E/Em , Em = 2eeg7 is the fractional neutrino energy. 


Dependence of the total rate on parameters eec /,7 taken as independent is oc Along with the 

relation Em = 2 eeg 7 , conclude that it is desirable to choose highly stripped heavy ions in order to achieve 
both high energy neutrino and large production rates. Rates further depend on N\peg{0)\'^ of injected ion, 
which requires a coherence of large peg{0)- We shall derive a constraint on this coherence factor in Section 
6 . A numerical estimate then gives 


F = V Fi ~ 3.1 X 102^Hz(-^)^/2 


S!N\p,giO)\ 

108 


^eg 


^50keV 




with Em = 2 eeg 7 


IGeV 


^eg 7 
50keV 104 ■ 


(37) 

(38) 


We may offer an interpretation of dependence of the total rate on involved various quantities. Ignoring 
dimensionless numerical values one has the relation, 


F oc i • N\p,g{0)\^ • GlEi -^g. (39) 

Each factor written here has a clear meaning. What this dependence implies is a scaling law with the 
boosted factor 7 in the laboratory frame of the circular motion. Note first that there is a hidden 7 factor 
in the radius of l/p = QeB/{'^MA) with Ma the ion mass and Qe the charge of ion. Except for the first 
factor 1/7 which arises from the prolonged lifetime oc 7 , other factors are dictated by the simple scaling of 
the basic atomic energy, with eeg l^-eg- This scaling law holds in the photon emission rate discussed in 
Section 6 . 

For a variety of expected neutrino experiments based on CP-even neutrino beam, it is important to have 
a beam of neutrino energy high enough in the GeV region (at minimum, larger than 0(200) MeV), since 
only then one can clearly detect the charged current (CC) interaction of If this requirement is not 

fulfilled, one only has CC interaction of and all kinds of neutral current (NC) interaction including 

NC process has a lower rate and experiments are harder. The GeV neutrino production requires 
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2eeg7 — IGeV for this combination of the ion parameter and the boost factor. For 7 < 10'^, it is necessary 
to have €eg > 50keV in order to reach 1 GeV neutrino energy. It is then important to excite electrons deeply 
bound in highly stripped ions in order to reach the keV binding energy of valence electrons in ion. This 
might be a non-trivial problem, but we assume that this is possible m- Judging from Fig([3]) it appears 
that there is an excellent chance of neutrino experiments in 0(0.5 ~ 1 GeV) energy range, which roughly 
gives a large rate of order, 10^° ~ lO^^Hz of the neutrino pair emission. A large value of eeg and a large 
boost factor 7 are required to accomplish this goal. 


Equi-rate, equi-energy curves 



keV 


Figure 3: Total equi-rate (given by ea. (1371) 1 curves for 3 x 10^^, 3 x 10^^ Hz’s in colored solid in (ce^/keV, 7 ) 
plane. Also plotted is Em = IGeV curve in dashed black and Em = 0.5GeV curve in dashed-dotted black. 
V|peg(0)P*S'g = 10® is assumed. The normal hierarchical mass pattern of the smallest neutrino mass 0 is 
taken, with egg = 50keV,7 = 10^ the synchrotron radius, p = 4km. 


Very importantly as a caveat, a new scheme of continuous injection or generation of coherent excited ion 
beam should be invented, because an excited ion, once it produces the neutrino pair, is not expected to be 
reusable for another source. The design and realization of this scheme might be challenging, and one may 
have to tolerate a sizable and unavoidable reduction of effectively usable ion number. 

A possible problem of highly stripped heavy ions is their large magnetic dipole (Ml) transition rate |14] . 
For concreteness let us take an example of He-like ion , Pb®*^"*". A good candidate for the initial ionic level is 
|e) = ((2s)(ls)) (a spin triplet state described in jj coupling scheme) of level spacing eeg ~ 70keV. The 
beam loss rate due to Ml photon emission is 

= 7MlA^/Oee(0) , (40) 

where the Ml decay rate 7 mi is ~ 3.4 x 10^®Hz according to |15] . By requiring that this loss rate is smaller 
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than the neutrino pair emission rate, one derives 

Ift5(0)|2 > O(0.1)p„{0)(^)“‘ ■ (41) 

Assuming the relation \pf,g{fi)\^ = /?ee(0)Pgg(0) that holds for a quantum mixture of pure states, one may 
further derive a lower bound for Pgg{0) of order 0 . 1 ( 7 / 10 '^)“^. This condition in the general case of He-like 
ions, with the inequality Pgg{0) < 1, further gives constraint on a relation of the boost factor and the level 
spacing. Details on these constraints should be worked out after more detailed R and D investigation on 
candidate heavy ions. 


6 High energy gamma ray beam 


For completeness we shall present main results for high energy gamma emission from the circulating excited 
ion. The high energy gamma ray emission occurs between different parity states among which El transition 
is allowed. The formalism in the main text is readily adapted to this case and we shall be brief in presenting 
results. 

The basic hamiltonian operator of El photon emission is (using a similar notation as in the text) 


= f d^x—A{x) ■ J^{x), J^{x) = ^7 j dtpeg{t)pegd{x - XA{t )), 

J y/l J 


(42) 


where A{x) = e^^'^e^ly/2ojV is the vector potential of emitted plane-wave photon (e^ being the polarization 
of photon). The other contribution arising from the center of mass (CM) motion part of ion as a whole (x Pa 
(atomic momentum) has been omitted, because it does not contribute to the internal atomic transition of 
1^) ^ \g)- The CM part gives contribution similar to the usual electron’s synchrotron radiation and gives 
rates much smaller than the rest of contribution. See more on this point. 

The matrix element of the internal part oc peg leading to ea. (l42p has been worked out as follows. The 
relativistic form of interaction hamiltonian density after the Lorentz boost of 7 factor is given by 

67 / d^x{g\(^A-a.'ijj + A-, (43) 


using the radiation gauge in the atomic rest frame, where a is the Dirac 4x4 matrix, a = 707 , and /3 
is the Lorentz boost vector. The orthogonality of (non-relativistic) wave functions of |e) and \g) gives the 
vanishing second contribution cc ijj to the leading first order to v/c {v being the velocity of atomic electron 
in its rest frame). The first contribution gives the internal contribution of eq. (j42p when the matrix element 
f d^x{g\e^^'^ip^aip\e) is written in the atomic rest frame, taking the long wavelength approximation of fc —>■ 0 
valid comparing with a larger inverse atomic length scale. The atomic matrix element may further be recast 
into the usual dipole form, using the equation of motion; Peglme = —ieegVeg with eveg the dipole matrix 
element for El transition. 

Calculation of the phase integral involves the energy-momentum (w, fc), |/c| = a; of a single photon. 
Stationary phase points appear due to the presence of the energy eeg in the phase integral. Straightforward 
calculations using the same approximation as in the text lead to photon emission rates. It would be 
instructive to start from a detailed discussion of the angular distribution. The double differential emission 
rate is given, to the best accuracy we know of, by 


d^T 

dwdVL 


(cosecos^)^^ -a,+ 1(1 - (44) 
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where the squared dipole moment f'^g was replaced by the related decay rate 'jeg (Einstein’s A-coefficient). 
The bracketed quantity in the argument of the negative fractional power —1/4 must be positive definite. To 
the leading order of the boost factor 7 , this quantity is approximated near the forward direction by 

7 272 2 ^ 27 ^ ^ 

The positivity requires the inside region of an ellipsoid in (0, ip) plane, 


^2 ^ £e£^2 < UJ 


'yuj 


72 ^ 


(jJm = 276 


eg 1 


(46) 


along with uj < ojm- Thus, there exists an interesting angular asymmetry: the cylindrical symmetry around 
the tangential direction is broken. With this approximation, the double differential rate becomes 


(fV 

diodn 


1 

2 ^/^ • IGtt^ 


A^|Pe,(0)p7VP— 







(47) 


The small angle approximation here is valid only for a small value of — ca)/( 7 ^a;). The approximation 
clearly breaks down at the infrared limit w —>■ 0 . 

Further angular integration is straightforward if one uses the small angle approximation, leading to the 
photon energy spectrum and finally the total emission rate. 


(£_ 

du 

r = 


;^iV|pe3(0)p7^^/^(^)^^^(l - ■ 

^eg 

12^2'^l^eg(0)|^7^7esV^, / = ^ dyy^/^{l-y)^/^ 


247r2' 

I 


~ 0.4165. 


(48) 

(49) 


The dimensionless spectrum function — x )3/4/, 

,x = ojjujm is plotted in FigQ after renormalization. 
Its end point is at ujm = ‘2'yeeg and the averaged energy value is 0.42wm- A typical value of the total photon 
emission rate is 



Figure 4: Normalized universal energy spectrum of photons emitted from excited ions. 


12 


















In a similar fashion as for the neutrino pair emission, the formula of the total El photon emission rate 
may be interpreted using the 7 —scaling law. Let us ignore dimensionless numerical factors for this purpose. 
The total rate then has dependence on various quantities: 

T (yi^-N\peg{^)\^ ■ e^reg'^Lol^-Um = ‘iieeg. (51) 

When one regards the atomic dipole ef^g as an invariant and intrinsic quantity to atom, other energy factor 
scales as oc 7 under the Lorentz transformation, along with the prolonged lifetime factor I /7 in front. This 
law explains 7 and e^g dependence of photon emission rate oc ^"^ejg ^ as well as that of neutrino pair emission 
oc {p regarded as 7 independent). 

Comparison with the usual synchrotron radiation may be instructive and also interesting. We can work 
out photon emission caused by ion circular motion in the ground state (or kept in the excited state) or simply 
an electron’s circular motion, using the same calculation technique as above. The basic hamiltonian arises 
from the omitted oc Pa term and the calculation is purely classical unlike the semi-classical approximation 
in the case of photon emission from an excited level. There is no lifetime related factor oc I /7 in this case, 
because the synchrotron emission is not a decay process. The result differs in an essential way from the case 
of excited ion, in that there is no stationary point of time integral. The phase integral in this case takes the 
form, 


dxh{x) cos + ^x) ^ . (52) 

There is no phase cancellation unlike in the case of excited ion. Instead, the exponential cutoff emerges for 
large 

We can finally derive in the large radius (p) limit a compact result for the energy spectrum and the total 
rate: 



dio 


= N 



2tt Q^a 1 


3 dvr 7^ 


/ 

I bjjuic Vs 




UJ 

4km N 7 


Wc = — 7 ^ ~ 

2 p 


, 7 ,n4km 
75eV ^)3- 

p 


r = ^ 8.4 X lo^^Hzg^:-^ 

y 2 8 p ’ ^ p 1019 104 


(53) 

(54) 


where Qe is the charge of ion. The value of the total rate given here corresponds to 1C ions equivalent to 
~ IO 19 ion numbers. The value p in the upper bound of ^ integral is estimated around 5 from the available 
angular area of 47r. This result is in a fair agreement with the standard results given in textbooks such as 
[9], considering the crudeness of matrix element estimate given here. We show the spectrum in Fig([5]) for 
the reader’s reference. 

It is interpreted that for electron synchrotron radiation the Zeeman splitting energy eB/{pme) = l/p is 
extended by the 7 ^ factor. Use of the internal atomic energy in our problem has two important effects: (1) 
larger energy spacing than the Zeeman splitting, and (2) kinematical power law cutoff at 2eegP rather than 
the exponential cutoff jme = 7^/p ~ 27 ^neV( 100 m// 3 ) for synchrotron radiation. 

One may work out a requirement on the coherence Peg{0) by demanding that the synchrotron radiation 
is not an obstacle against the neutrino pair emission m- It is imposed that the number of emitted neutrino 
pair per revolution of circular motion T x 2 'nplc (equivalent to the number of de-excited ion caused by 
neutrino pair emission) is much larger than the number of emitted synchrotron photon per revolution. This 
condition gives a constraint on the coherence, 


Ift,(0)i»1X io-‘(3(m)-3/=(£-)-=/-‘{ 


ceg 


50keV' 


n-11/4 


(55) 


taking the spin factor to be unity. Si = \. If this condition is violated in the case of a large Q, one may 
have to think of compensating the loss of excited ions by irradiation of laser each time of revolution m- 
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number, power 


Synchrotron radiation spectrum 



Figure 5: Normalized synchrotron spectrum shapes. Number spectrum in solid black and power spectrum 
in dashed red. 


We have several comments based on results of the gamma emission from excited ions. First, a different 
kind of coherence effect over a larger volume may further enhance photon emission rates by the super¬ 
radiance mechanism of Dicke m- In the case of two-photon emission mentioned above the macro-coherent 
paired super-radiance (PSR) may further enlarge the coherent region |11] not restricted to an area of the 
photon wavelength in the Dicke case. It might even be possible to produce coherent gamma ray ’laser’, with 
a help of macro-coherence. As an example, 2s —Is two-photon transition of H~ ion may be an excellent 
source of coherent two-photon emission due to its long lifetime of 2s excited ion. The achievable energy is not 
large, however, of order 200 keV y/lO^ for hydrogen ion. Molecular vibrational transition u = 1 —)• 0 of pH^ 
may be better due to their easiness of Raman excitation. Recently, the macro-coherent PSR of vibrational 
transition of neutral pH 2 was observed [l 8 j, in which we achieved a macro-coherence of ~ several % over a 
target of 15 cm long. Rates of two-photon emission from circulating excited ions may be worked ont as in 
the rate calculation of neutrino pair emission. Our 7 —scaling law suggests that two-photon emission rates 
are large despite their effective, weaker coupling. 

Even as a technical strategy towards high intensity neutrino beam, it would be wise to first study basic 
experimental feasibility of heavy ion excitation aiming at high energy photon emission, since it would be 
easier to detect and study the mechanism of photon emission in detail. 

In summary, a new method of prodncing CP-even coherent neutrino beam from circulating excited ions 
was proposed. Large production rates of neutrino energies extending to much beyond the keV region were 
derived. When ions are excited to a different El allowed level, they may provide high intensity gamma ray 
beam much beyond the keV range. Evidently, much R and D works, both theoretical and experimental, are 
needed to determine a realistic design using a specific ion. 
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